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Abstract 

This paper investigates the performance of the downlink channel in distributed antenna systems. 
We first establish the ergodic capacity of distributed antennas, under different channel side information 
(CSI) assumptions. We consider a generalized distributed antenna system with N distributed ports, 
each of which is equipped with an array of L transmit antennas and constrained by a fixed transmit 
power. For this system we calculate the downlink capacity to a single antenna receiver, under different 
assumptions about the availability of the channel states at the transmitter. Having established this 
information theoretic analysis of the ergodic capacity of distributed antenna systems, this paper also 
investigates the effect of antenna placement on the performance of such systems. In particular, we 
investigate the optimal placement of the transmit antennas in distributed antenna systems. We present 
a fairly general framework for this optimization with no constraint on the location of the antennas. 
Based on stochastic approximation theory, we adopt a formulation that is suitable for node placement 
optimization in various wireless network scenarios. We show that optimal placement of antennas inside 
the coverage region can significantly improve the power efficiency of wireless networks. 

I. Introduction 

The main challenge in the design of future cellular wireless networks is the mitigation of the 
adverse effects of interference. One possible strategy to alleviate interference, both in the uplink 
and the downlink of cellular networks, is to reduce the overall transmit power by using distributed 
antenna systems (DAS), which also have the additional advantage of improving capacity and 
coverage [1]. Moreover, by reducing the access distance between the transmitter and the receiver, 
distributed antenna systems have direct impact on the energy efficiency of cellular networks, 
which may lead to greener architectures in the future. 



1 



Downlink capacity of distributed antenna systems can be studied in the context of the MISO 
wireless channel with a power constraint on each distributed antenna port. Ignoring this per 
antenna power constraint for the moment, the capacity of wireless MISO channels has been 
studied extensively in the literature, under different assumptions on availability of the CSI at 
the transmitter [2], [3]. For the case of having CSI both at the transmitter and the receiver, we 
know that beamforming is optimal [4], whereas for Rayleigh fading channels when the CSI is 
only available to the receiver, the optimal strategy is to send independent signals with equal 
power [2]. The main assumption in these studies is that we have a sum-power constraint across 
all transmit antennas. However, in a distributed antenna system, the transmitter cannot allocate 
power arbitrarily across the antenna ports, as each port has its own power budget that can be 
allocated only to its own antenna array. As a result of this per-port power constraint, the simple 
eigenvalue decomposition analysis of the MISO channels in [2] cannot be applied in this case 
to obtain the closed form solution. However, many numerical solutions have been proposed for 
dealing with restricted power constraints in the context of MIMO/MISO wireless channels. The 
downlink capacity of multiuser MIMO channels with a per antenna power constraint has been 
investigated in [5]. Also, based on the duality of the uplink and downlink, an iterative algorithm 
has also been proposed in [6] for solving the sum rate maximization problem by utilizing the 
convexity of the dual problem. None of these algorithms, however, provide a closed-form solution 
to the capacity, because of the complexity of the optimization problem. 

In this paper, we establish a closed-form solution for the capacity of the downlink channel 
for a distributed antenna system. We assume there is a single user that receives the signal from 
multiple transmit ports, each of which is equipped with an antenna array and a dedicated power 
budget. The special case of this problem has been studied in [7], where each transmit antenna 
port has a single antenna. We generalize the results in [7] for the case where each antenna port 
has multiple antennas and the antenna ports are geographically distributed in the cell coverage 
region. We consider two cases for availability of the CSI. In the first scenario, we assume that 
CSI is available both to the transmitters and the receiver and we show that the optimal signaling 
for each of the antenna ports is beamforming. In the second case, we assume that we have a 
symmetric fading channel, and the CSI is only available to the receiver. For this case we show 
that each of the antennas of each port should send independent signals. In other words, we show 
that the possibility for cooperation among the antenna ports does not help in increasing the 
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capacity, when we have a symmetric fading channel. For both cases, by solving the associated 
stochastic optimization problem, we establish the capacity in closed form. We also determine 
the impact upon capacity of interference from other cells, assuming this interference is treated 
as noise. 

Having established this information theoretic analysis of the ergodic capacity of distributed 
antenna systems, in the second part of this paper we investigate the effect of antenna placement on 
the performance of such systems. The main reason for looking into the placement problem is the 
fact that the performance gain of distributed antenna systems (the capacity gain or equivalently 
the power saving gain) is largely influenced by antenna locations. In [8], it was shown that 
the placement of distributed antennas can significantly influence system performance and that 
this optimal placement magnifies the advantages of using DAS over the traditional systems 
with centrally-located antennas [9]. In other words, the capacity increase of distributed antenna 
systems is largely influenced by antenna locations. However, many of the studies on the placement 
optimization of antennas [10]— [12] impose restrictions on the topology of the network, such as 
linear cells or antennas deployed along a circle. To address the problem without imposing these 
topological constraints, this paper investigates the optimal location of the transmit antennas in 
distributed antenna systems in a general framework. While our framework can incorporate a 
broad set of performance metrics, our results will focus on the capacity and power efficiency 
gains obtained through optimal placement. 

The remainder of this paper is organized as follows. First in Section II, we provide a basic 
description of our distributed antenna system. In Section III, we establish the capacity results for 
both cases of having channel side information at the transmitter or not. We also analyze the effect 
of interference from the neighboring cells on our capacity results in Section III-B. In Section IV, 
after formulating our placement optimization problem, we introduce two different approaches 
for finding the optimal placement. We investigate some power saving techniques that can be 
applied to our generalized DAS system in Section V. Section VI illustrates the performance of 
the proposed algorithms in different settings. Finally, in Section VII we conclude the paper and 
present some possible extensions. 
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II. System Model 

In this section we describe the system model for analyzing the downlink performance of 
distributed antenna systems in a cellular setting. Figure 1 is the general architecture of the 
distributed antenna system, DAS(N,L), in the multicell environment we consider in this paper. 
We assume that each cell is covered with a total number of N distributed antenna ports and 
each port has L microdiversity antennas. We assume that the antenna ports are all connected to 
the central antenna port via dedicated links (e.g fiber) with unlimited bandwidth and zero delay. 
We consider a one-tier cellular structure in which a given cell is surrounded by a continuous 
tier of six cells that can cause interference. The central cell in Figure 1 is indexed by j = 0, 
while the surrounding cells are indexed by j = 1, 6. We denote the positions of the ports in 
the central cell with a 1 x 2N vector P which is defined as 

P=[Pi,p 2 , ,Pat], (1) 
where each of the p/s indicate the position of the i th port in the Cartesian plane as Pj d = [x^ yi\ 

■ def 

for i e {1, 2, N}. Denoting the centroid of the jth cell in the Cartesian plane with o 1 = 
[o£, o£], we assume that the relative positions of the ports in the jth cell with respect to a? is 
identical for all j = 1,2,. ..,6. In other words, if we denote the location of the i th port in the 
jth cell with p^, we assume that for all antenna ports (i = 1, 2, N) in the neighboring cells 

( j = 1, 2, 6) we have 1 

Pl-Pi + O 7 (2) 

All of the NL antennas in our N distributed antenna ports together construct a macroscopic 
MISO system with the NL x 1 channel vector 

Kj) = [hi(3)M3),-,hN{j)l (3) 

where h n (j) denotes the channel vector gain from the L antennas of port n to the user at location 
u which is a function of the distance r between the nth transmitter port and the user. Because 
of the different distances between the antenna ports and the user, our channel vector gain should 
encompass not only the small scale fading but also the large scale fading together with the path 

'Note that we put the origin of the Cartesian plane at the centroid of the central cell (j = 0) 
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loss as 



K{j) 



A3) 
iJn 



Aj) Ai) Aj) 

Jn,l' Jn,2i J n,L 



(4) 



where f^J denotes short term fading in the channel at the ith antenna of port n. In particular, 
we assume that the fading coefficients are symmetric in the sense that the distribution of —f^\ 
is identical to the distribution of f^J for all 1 < n < N, 1 < i < L. One example of such 
a symmetric fading channel that can be readily applied to our formulation is Rayleigh fading, 
where f^J is modeled by a complex Gaussian random variable distributed as CA/"(0, 1). Also 
g$ is a random variable representing the shadow fading. Cn\r) in (4) is the power path loss 
function which is a function of distance r and also depends on propagation frequency, as we 
will introduce later. In this paper we consider an exponential path loss function of the form 
C(r) = (5 x r a , where a is the path loss exponent whose value is normally in the range of 
two to six. We assume independent shadow fading between each port and also independent fast 
fading at each antenna per port. Here, the distance function r(-) can be calculated as 

r( Pi ,u) = max{£>( Pi ,u),r 0 } i e {1, 2, N}, (5) 

where ©(Pj, u) is the Cartesian distance between the user at location u G II and the i th antenna 
port's location Pj and can be written as 

^(P;, U) = \j{P X -U X ) 2 + {Py-Uy) 2 . (6) 

Also, r 0 in (5) is the minimum allowable value of di for which the far field approximation of 
the antenna propagation model is valid. 

Let the transmit signal vector of the jth cell be x(j) = [xi(j), x 2 (j)---, xnU)], m which x n (j) 
denotes the 1 x L transmission vector of the nth port of the jth cell. Then the received signal 
of the user in the central cell can be written as 

6 

y = h(0) T x(0) + J2V^J h U) T x(3) +n,= h(0) T x(0) + z, (7) 

where n is an additive Gaussian noise with variance a 2 and 7j is a positive real multiuser coding 
gain that quantifies the fact that the signals from the jth cell are not completely orthogonal to the 
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signals of the central cell 2 . In deriving (7), since the neighboring cells are sending independent 
messages, we approximated the sum of interfering signals plus noise as a complex Gaussian 
random variable z with variance <j 2 z . We also assume that each port has a separate transmit 
power budget of S n (n — 1, 2, N), i.e 



E 



Xn(j)x n (j) H = S n , (8) 



for all ports 1 < n < N, and in each of the cells 0 < j < 6. Note that in (8), the nth transmission 
port can allocate its own power budget S n among it's antennas, but different ports have no ability 
for allocating power amongst themselves. We also denote the transmission power constraint of 

— def 

all ports by the vector S = [S±,S2, ...,Sn]- 

III. Downlink Capacity of Distributed Antenna System 
A. On downlink ergodic capacity 

Given the channel model in Section II, for fixed port location matrix P and fixed user's location 
u, we can write the average ergodic rate of a single user under different circumstances. In this 
section, we consider two cases of channel information at the transmitter: the channel coefficient 
vector h(j) is known to the transmitters of the jth port (channel with CSIT) or not (channel 
with no CSIT). In both cases, we assume that the channel coefficient vector is known at the 
receiver, which is a good approximation for practical systems with receiver channel estimation. 
The capacity of this channel depends on the power constraint on the input signal vector x in (7) 
and also on the availability of the CSI at the transmitter. The goal of this section is to investigate 
the downlink average ergodic capacity of a user in the following scenarios for different sets of 
information that are available at the transmission ports. 

1) CSI only at the receiver: In this case, we assume that the channel vector h is unknown to 
the transmitter. Because of the Gaussian noise and known channel at the receiver, the optimal 
input signal is Gaussian with zero mean [2], [7], and the achievable average transmission rate 

2 This coding gain can be set to one if we do not use orthogonal codes in the neighboring cells. On the other hand, = 0 
indicates that we have no interference from the neighboring cells. Also in a CDMA system, this can quantify the spreading 
code gain. 
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of the user for this case can be found by solving the following stochastic optimization problem 

1 



Max : E f 



log 2 1 



at 



h(0) Q h(0)") 



(9) 



Subject to : Q y 0 , tr (q nn ) < S n , VI < n < N , 

where Ef [•] denotes the expectation with respect to both slow and fast fading, Q = E x(0) H x(0) 
is the NL x NL Hermitian covariance matrix of the Gaussian input of the central cell, and q mn 
is the part of the covariance matrix Q that corresponds to the correlation of the transmitted 



x m (0) T x n (0) 



Note that the constraint 



LxL 



signal of the mth and the nth ports, i.e. q n 
tr(?Tm) < Sn corresponds to the power constraint of each antenna port defined in (8). The 
main issue in solving (9) is finding the optimal covariance matrix Q. Since the per-port power 
constraint tr (q nn ) < S n is not equivalent to a constraints on eigenvalues of Q, the analysis of 
[2] cannot be applied here. In other words, after the eigenvalue decomposition of Q — U A U T , 
although h(0) U has the same distribution as h(0) , the partial constraints on the eigenvalues of A 
would not be the same as the constraint on (9). Therefore, the problem is no longer equivalent to 
the analysis of the fading channels in [2] through eigenvalue decomposition. However, inspired 
by the work in [7] and based on the symmetric distribution of the fading channel, we can show 
that the optimal solution of (9) is Q* csir = diag(^lL, ^l, l^-l), where l L denotes anLxI 
identity matrix. This optimal solution means that each antenna port should transmit independent 
signals at full power. Somewhat surprisingly, our optimal transmission strategy in this case is 
similar to the transmission strategy of a multiple access channel. In other words, for a symmetric 
fading channel, e.g. a Rayleigh fading channel, the possibility for cooperation among the antenna 
ports does not help in increasing the capacity. In this case the capacity is 

1 



C(u, P, S) = E f 



log 2 1 



M0)Q:. Si M0) H ) 



E 



/ 



N 



z n=l 



(10) 



where |.| denotes the / 2 -Norm. 

2) CSI both at the transmitter and the receiver: In this part, we assume that the channel 
vector gain h is known at the transmitter. Similar to the case of having CSIR, the optimal input 
signal is still Gaussian with zero mean, because of the Gaussian noise and known channel at 
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1 



Max 
Subject to 



log 2 1 + - h(0) Q h(0) H ) 



(ID 



Q ,tr(q nn ) < S n , Wl<n<N 



In [7], the problem is solved for the special case of L = 1 by relaxing the positive semi-definite 
constraint Q y 0. Here, we find the solution for a general DAS(N, L) system. We show in 
Appendix II that the optimal covariance matrix Q* csit is 

Qlsit = [Qi,Q2, Qn] H [qi, 92, Qn] , (12) 

where q n d = V~Sn, I < n < N. Also, the optimal signaling solution for the nth antenna 

port is beamforming with a beam weight of h* (0)/||h n (0)||. As we can see, in each antenna port, 
the power allocation among L antennas is proportional to the amplitude of the channels, while 
the phases are matched and the total transmit power is set to its maximum allowed budget. For 
the special case of L — 1, each antenna port should match the phase of its associated channel 
so that the signals add coherently at the receiver and also set the amplitude of the transmitted 
signal independent of the channel and fixed according to the power constraint. In other words, 
the key information that we need for achieving capacity is the phase of the channel gains, but 
not the amplitudes. Note that for the case of L = 1, there is actually no power allocation among 
the transmit antennas and the transmit power of the nth antenna port is fixed at S n . On the 
other hand, if we only have a single antenna port (N — 1) with L microdiversity antennas, 
beamforming with the weights that are proportional to the channel gains is essential to achieve 
the maximum capacity. 

The average ergodic rate of DAS(N,L) system with CSIT can then be written as 

1 



C(u, P, S) = E f 



log 2 l + -h(0) Q* csit MO)") 



E, 



z 
N 



. n=l 



(13) 



where Ef [•] denotes the expectation with respect to both slow and fast fading. Comparing (13) 
and (10), we see that under the per-port power constraint, the presence or lack of channel phase 
information at the transmitter has a significant impact on the optimal transmit strategy and hence 
on the channel capacity. 
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B. Interference from neighboring cells 

We also investigate the effect of interference from the neighboring cells on ergodic capacity 
of the system both with and without CSI. We assume there is no cooperation between cells 
and that intercell interference is treated as noise. As a result, regardless of the availability of 
CSIT in a given cell, the interfering signals from neighboring cells add incoherently at each 
user's receiver. Therefore, we can write the power of interference plus noise (complex Gaussian 
random variable z in (7)) for a user in the central cell (j = 0), as 



where we assume that all antenna ports in all neighboring cells are using the same transmit 
power of S. Therefore, by replacing a 2 z with the variance of the noise in (10) or (13), we 
can consider the effect of interference on the expected ergodic capacity. 

In order to illustrates the effect of interference on the downlink capacity of DAS (7, 1), we 
consider the following model. As shown in Figure 1, we assume that one port is in the center 
and six other ports are in a circle of radius r = R/2 around it, where R is the radius of the 
hexagonal coverage region. We assume independent log-normal shadowing of variance a sh = 8 
between each port and the user. The path loss exponent is assumed to be a = 4 and we set 
the power constraint on each port 3 such that the SNR of a user at distance R is 10 dB. Figure 
2 illustrates the downlink average ergodic rate of a user, averaged over both the shadowing 
distribution and user location which is assumed to be uniformly distributed inside the cell. We 
consider both cases where CSI is only available to the receiver and where CSI is available both 
to the transmitter and the receiver. As expected, having the CSI at the transmitter improves the 
average ergodic rate substantially. It is also worth mentioning that in the case of having only 
CSIR, the SINR would be much worse for the users on the edge of the cell in comparison with 
the case where we have CSIT (since the interference terms will add incoherently even when we 
have CSIT, whereas the signal is adding coherently). 



In this section we focus on the placement optimization of the antenna ports in DAS(N, L). 
As we will show later in Section VI, depending on different assumptions about the availability 

3 in absence of any interference from other ports 
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N 




(14) 



IV. Antenna Placement 
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of CSI, the optimal antenna placement can have a huge impact on the energy efficiency of 
distributed antenna systems. For the sake of simplicity, we neglect the effect of fast fading in 
our optimization problem, i.e we replace the fast fading terms in the channel gain vector (4) with 
f nji — 1, for all 1 < n < N, 1 < % < L. In addition, we start with the scenario when there is no 
interference and the metric that we want to optimize is the cell averaged ergodic capacity. The 
averaging here is done with respect to the location of the user, which is uniformly distributed 
within the cell region. We use u to denote the position of the user. The ergodic capacity in 
bit/ s/ Hz conditioned on the user's position u G II, antenna placement P 6 II and also the 
power budget S is introduced in Section III for two different cases where the CSIT is available 
at the transmitter or not. Specifically, depending on the use of (10) or (13) for the capacity, the 
cell averaged ergodic capacity C for a particular placement of the antennas can be calculated 
by averaging C(u, P, S) over user positions u inside the cell region II as 

C(P,S)=E U C(u, P,S), (15) 

where E u {-} denotes the expectation with respect to the location of the user and in our model 
we assume that it has a uniform distribution inside the cell region. 

Now we can describe the placement optimization problem for DAS that we consider in this 
section. 

Problem 4.1: Consider N to be a positive integer quantifying the number of ports in our 
distributed antenna system. Let the variance of the log-normal shadowing a sh and the minimum 
distance for far-field approximation r 0 be given design parameters. Also assume that the power 
path loss function £(•) is known and we have the same power budget at each antenna port, 
denoted by S, i.e S = [S,S, ...,S]. We want to find the solution to the following optimization 
paradigm: 

P* = argrnin,S, (16) 

subject to : 

C(p,S)>C t , (17) 

p t eU ie {1,2,..., N}, (18) 

where C t is the target for the expected cell average ergodic rate. Here P, as defined in Section 
II, is the matrix containing the collection of locations for all the ports in the system. 
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Remark 4.1: It can be easily seen that Problem 4.1 is equivalent to maximizing the expected 
cell average ergodic rate with a constraint on the maximum power budget S. 

In the remainder of this section, we propose two different solutions for the aforementioned 
placement optimization problem. First we present a simple approximation method that casts the 
problem in a simpler space that is independent of channel characteristics. Next we propose a 
very general framework for placement optimization that can handle more complex optimization 
metrics. 

A. Simple approximation method 

Either of the ergodic rate formulas (10) or (13) introduced in Section III involve a sum over NL 
nonnegative variables among which the L terms of the antenna port with the smallest distance 
r dominates the sum. One approach for tackling the placement optimization problem is to only 
consider the dominant terms in the summation and ignore the rest. Starting with (10) where the 
coherent transmission is not possible due to the lack of phase information at the transmitter, by 
considering the closest antenna port to the user we can approximate the average capacity with 
the following lower bound [10] 

C(P,S)>E U E f \og 2 (l + -^-) : (19) 

where 

r min = min V(Pi,u), (20) 

l<i<N 

m = arg min X>(p,-,u). (21) 

l<i<N v 1 ' 

When CSIT is available, the approximations in (19) remain unchanged, except the fact that the 
power term S should be multiplied by L, because of the coherent summation of the signals at the 
receiver. However, this constant L has no effect on our placement optimization problem so we 
continue with (19) as a lower bound on average capacity. Also, note that the approximation in 
(19) works better for higher values of path loss exponent a, since the larger values of a magnify 
the difference between the dominant term and the second largest value in the summations of 
(10) and (13). By exchanging the order of the two expectation operators and noting that the 
function f(x) = log 2 (l + |) (a > 0) is a convex function, we can apply Jensen's inequality to 
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obtain a lower bound for averaged capacity as follows 

C,P, S) > E, E. ,og 2 (l + J^)> E, .o g2 (l + .J^) . (22, 
We can use a similar approximation, based on Jensen's inequality, to average over shadowing. 

def 

To this end, we define the auxiliary random variable g m = log 10 g m and rewrite (22) as 

= l og2 (l + — | -Y (25) 

Recall that, according to our model, the shadowing term g m has a log-normal distribution. Hence, 
g m in (24) has a zero mean normal distribution. Also, the inequality in (24) is a direct application 
of Jensen'e inequality to (23) by noting that the function log(l + a 1CP) is a convex function 
of x e R, for any positive coefficient a. Note that equations (23), (24) and (25) suggest that 
we can ignore the effect of shadowing and still obtain an upper bound on the minimum power 
required to achieve a certain average ergodic rate of the cell. We will confirm this observation 
in Section VI via simulation, in which we illustrate that shadowing has a negligible effect on 
the placement optimization of the ports. 

Now, in order to obtain an upper bound for the solution of Problem 4.1, we can replace 
constraint (17) with a more conservative constraint of the type 

log 2 (l + = f/ .f u ) > C t . (26) 

V E u {C(r min )} J 

This inequality indicates that the approximate solution 4 of Problem 4.1 can be obtained by 
minimizing E u {C(r min )} . Therefore, with this lower bounding we can get a criterion for antenna 
location design which is to minimize the expectation of the path loss function C (for any such 
function) from the randomly distributed user to the nearest antenna port. The following example 
illustrates this bounding approach for a simple example, where we want to optimally place three 
ports within a circular coverage region. 

Example 4.1: Let r 0 be zero, the number of ports N = 3 and the cell region be a hexagon of 
radius R, as illustrated in Figure 3. The goal is to find the approximate solution of Problem 4.1 

4 note that replacing (17) with (26) will shrink the feasible set of Problem 4.1 and hence give an upper bound on the solution 
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for this setting. Due to the symmetry of the problem, one can argue that the optimal placement 
of these three single- antenna ports should also be symmetric around the center of the region. One 
possibility for such a symmetric placement is illustrated in Figure 3 where the regions separated 
by dashed lines are sectors with central angles of 2tv/3 and the antenna ports are located on the 
bisector of the central angle. Here, for a randomly placed user, r min in (20) would be the distance 
between this user and the associated port in its sector. Now, in order to find the optimal location 
of the antenna ports, we need to minimize E u {£(r min )} or equivalently E u {r° in }, where the 
expectation is over a randomly placed user u within the sector. The solution is the centroid of 
the sector. The centroid here depends on the path loss exponent a (for instance, for a = 2 this 
reduces to the geometric centroid). Using this lower bound approximation, the optimal radius r* 
to deploy the antenna ports ranges from r* = 0.57 R to r* = 0.58 R, depending on the path loss 
exponent a (see Figure 5). This indicates that, using this lower bound approximation approach, 
the optimal position does not change much for different values of a. Note that this solution is 
based on the aforementioned approximation approach and has some implicit assumptions and 
approximations that make it suboptimal. 

As we can see, the antenna port locations obtained with this approach only depend on the 
antenna port number N, the user distribution (including cell shape), and the path loss exponent 
a. However, these optimal locations are completely independent of the availability of the CSI 
at the transmitter. In other words, using either (10) or (13) for the ergodic rate yields the same 
solution. Moreover, the solution is also not affected by the underlying shadowing model in the 
system. Note that the proposed criterion aims to maximize a lower bound of the performance 
metric C, rather than C itself. Therefore, the obtained antenna port locations are suboptimal. 
It is also not possible to generalize this model in order to take into account the interference 
from other users within the cell or from neighboring clusters. The other problem is that this 
approach has an implicit assumption that in the optimal placement, all the ports would be placed 
far away from each other so that the lower bound (19) is tight, which might not be true in 
general. However, this method becomes more accurate for the regime where we have high path 
loss since the lower bound (19) becomes tighter. 

In the next subsection, we propose a more general technique based upon stochastic approx- 
imation for optimizing the placement, that can handle a much more general class of objective 
functions. 
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B. Stochastic Update of the Locations 

In this part, based on stochastic approximation theory, we introduce a general framework for 
solving the placement optimization paradigm. For the sake of clarity, rather than solving (16) 
directly, we will first focus on the equivalent problem of maximizing the average ergodic rate 
with the given power budget S. We touch upon the generalization of the proposed algorithm for 
more sophisticated optimization metrics and also the direct solution to (16) afterwards. 

Assume we want to find the optimal placement vector P in order to maximize 

C(P,S)=E U E f C(u,P,g,S), (27) 

def 

where g = [gi, g^ <?jv] is defined to be a random vector of size N, indicating the shadowing 
and E/[-] denotes the expectation with respect to this shadowing. Also C(u,P, g, S) is the 
instantaneous ergodic rate for a given realization of the user position u, shadowing vector g, 
and given location matrix P, depending on which information we assume to be available at the 
transmitter. In particular, when we assume that CSI is only available to the receiver but not the 
transmitter, from (10) we can write the capacity as 

and when CSI is assumed to be known both at the transmitter and the receiver as in (13), the 
instantaneous capacity can be written as 




Now, in order to find the optimal location matrix P* that maximizes (27), we use the Robbins- 
Monro procedure [13] from stochastic approximation theory. Before getting into the algorithm 
that finds P*, we will briefly describe the Robbins-Monro procedure. In stochastic approximation, 
maximizing a utility function U (x) can be done via an iterative update of the optimization variable 
x in the following way 

x t+1 = x t + a* (U x (x t ) + M t ) , (30) 

where U x (x) d = M t is some zero mean noise in the estimate of U x (x) and {cr*}^ 1 is a 
square summable but not summable sequence. 5 In our specific placement optimization problem, 

5 The convergence of this Robbins-Monro iteration is shown under broad conditions in and [13, chapter 6]. In particular, the 
iterations converge to a local optimum when the noise terms M t have bounded variance. 
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we can look at the expected ergodic capacity as the utility function that we want to maximize 
and use the Robins-Monro procedure to update the location of the ports at each iteration. 
We can now specify the following iterative algorithm for solving (27). 

1) Initialize the location of the ports randomly inside the coverage region and set t — 1 

2) Generate one realization of the shadowing vector g t , 6 based on the probabilistic model 
that we have for the shadowing of the channel. 

3) Generate a random location u t for the position of the user based on the distribution of the 
users in the coverage area. 

4) Update the location vector as 

P t+ i = P t + <7 — | Pt (31) 

5) let t — t + 1 and go to step 2 until convergence. 

In order to provide some intuition about this algorithm, we will try to put it in the form 
of the Robins-Monro recursion (30). As we can see, dc( -»^s)) is an 

unbiased estimator of 

E u E f dc{u d /' g) . Hence by defining 

-0C(u,P,g,S) 



Mt ^cM,s) _ EuEf 



(32) 



<9P 

= dCiUt f p g - B) - ^{E u E f [C(u,P,g,5)]}, (33) 

we can rewrite the position update (31) of step 4 in the form of the Robins-Monro recursion by 
noticing that in our placement problem the utility function £/(•) = E u EfC(-) and the zero mean 
noise term M t is defined in (33). 

Now we describe how we can solve (16) directly, rather than solving the equivalent problem of 
maximizing the average ergodic rate with given power budget S. Similar to [14], a primal-dual 
algorithm can be used to solve (16) with the same stochastic update mechanism. Briefly, the dual 
problem and the Lagrange multiplier can be updated in an iterative fashion with a stochastic 
gradient method [14]. Unfortunately, due to the non-convexity of the objective function, there is 
no guarantee that the algorithm converges to the globally optimum solution, and it's possible that 
stochastic gradient steps may converge to a local optimum. This problem has been addressed 
for the special case where the path loss function has an exponential fall off with distance, where 

6 We can also consider the fast fading vector f t in our gradient estimation 
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we showed the placement optimization problem can be reformulated as a convex optimization 
problem [15]. 

In general, we do not have to restrict ourselves to the expected ergodic capacity metric 
E u EjC(-) in (27) and it can be replaced with other metrics such as outage probability. This 
can be done simply by replacing C(-) in (31) with any arbitrary performance metric that we 
want to optimize. In particular, in the following part, we show how the interference from the 
neighboring cells can be considered in the antenna port placement optimization. 

C. Interference Effect 

In this part, we investigate the effect of interference from the neighboring cells on antenna 
port placement optimization. As we mentioned in Section III-B, by replacing a 2 z from (14) with 
the variance of the noise a 2 , in (28) or (29), we can consider the effect of interference on the 
expected ergodic capacity. The same stochastic update algorithm of Section IV-B can be applied 
here for placement optimization of antenna ports, simply by replacing the capacity formula in 
(31). The only difference here is the assumption that we have about the location of the ports 
in the neighboring cells. As we described in Section II, we assume a similar relative position 
of all antenna ports in each cell with respect to the center of that cell. Therefore, by updating 
the location vector in (31), the location of all interfering ports would also change accordingly. 
This indeed forces the final optimal layout of the antennas (after convergence of the algorithm) 
to be identical in all the cells j = 0, 1, 6. As we will show in the next section, the effect of 
interference in the optimal location of the antenna ports is to move the ports towards the center 
of the cell in order to minimize the power of the antenna ports outside the cell boundaries. 

V. Power saving techniques 

In this section we study some power saving techniques that can be applied to our generalized 
DAS system. In particular, in the first part, we study the optimal output power of the antenna 
ports with respect to their position in the cell in order to mitigate the effect of interference. In the 
second part, we analyze the single transmission strategy where only one antenna port is used in 
the downlink. We also generalize the stochastic gradient method of Section IV-B in order to find 
the joint optimal location and power allocation of the antenna ports for the single transmission 
strategy. 
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A. Distribution of power 

Until now, we always assumed that we have the same power constraint S on the output 
power of antenna ports in the cell. The alternative is to adjust the power constraints of the ports 
according to their position in the cell. In other words, instead of having a fixed power vector 
S = [S,S, ...,S], we can treat the power budget of each antenna port S n (1 < n < N) as an 
optimization variable. In order to have a fair comparison with the case that all of the output 
powers are fixed to be S, we need to put a constraint on the total consumed power in the cell 
as J2iLi Si = N x S. Hence the averaged ergodic rate maximization problem becomes 

Maximize : C( P, S) (34) 

N 

Subject to : ^ S n < N x S, 

n=l 

where P and S are optimization variables. A simple variation of the stochastic gradient method 
of Section IV-B can handle this constrained optimization problem ( The algorithm is not included 
due to space limitations.) As we will observe in the numerical results of the next section, by 
increasing the effect of interference (increasing 7) the optimal power allocation algorithm tends 
to dedicate more power to the central antenna ports in order to reduce the amount of leaking 
power to the neighboring cells that cause interference. 

B. Single transmission strategy 

Several transmission strategies are possible in a distributed antenna system [16]. Until now, we 
assumed that all of the antenna ports are used in transmitting the signal, either in a coherent or 
non-coherent fashion. The alternative is to use the closest antenna port to the user for transmitting 
the signal. This way we minimize the required transmit power that also results in reducing the 
interference caused to other cells. As we will show in the numerical results, for a fixed average 
ergodic rate target, the single transmission mode increases the power efficiency of distributed 
antenna system significantly. 

VI. Numerical Results 

In this section we illustrate the performance of DAS and different algorithms we introduced 
in Sections IV and V. We start with the placement optimization algorithms of Section IV and 
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continue with different numerical studies for illustrating the power efficiency and area spectral 
efficiency of distributed antenna systems, under different assumptions on availability of CSI at 
the transmitter. 

A. Convergence of the algorithm 

In this part, we mainly focus on the convergence of the stochastic update method. We consider 
a distributed antenna system with N = 3 ports, each equipped with a single antenna. We assume 
that the coverage region n of these three ports is a hexagon of radius R = 1000. We assume 
independent log-normal shadowing of variance a sh = 8 between each port and the user. The path 
loss exponent is assumed to be a = 6. We initialized the location of the antenna ports randomly 
inside the cell coverage region and the initial points are marked with yellow circles in Figure 
4. Then we generated realizations of the random shadowing vector f t together with the random 
location u t uniformly inside the hexagonal cell region. We applied the algorithm in Section IV-B 
and plotted the trajectory of the locations of three ports as the algorithm evolved in Figure 4. The 
algorithm ran for 200k cycles. As we can see the optimum placement is symmetric with respect 
to the center of the cell region and they are all located in a hexagon of radius r opt = 550. This 
radius is checked to be independent of the initial placement of the ports by many repetitions of 
the algorithm with different initial values. 

B. Comparison of different optimization approaches 

In this part we compare the performance of the methods proposed in Sections IV-A and IV- 
B. Figure 5 illustrates the optimal radius r to deploy three ports, obtained from two different 
approaches introduced in Section IV. The plot illustrates the optimal placement for different 
values of path loss a. We used the stochastic update method for two different scenarios described 
in Section III and plotted the optimal deployment radius as a function of path loss. As we can 
see, depending on using (13) or (10) for capacity, the optimal placement strategy is different, 
especially for small path loss exponent a. In other words, different assumptions in Section III 
for the cooperation among the ports and also availability of CSIT at the transmitter will change 
the optimal placement strategy. Note that these assumptions do not affect the lower bound. As 
we can see, for larger values of a, the lower bound approximation approach creates a better 
estimate of the optimal deployment radius. This is understandable since the approximation we 
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used in our lower bound optimization approach is tighter in high path loss regime. On the other 
hand, even for low path loss exponents, the stochastic update method converges to the optimal 
solution. In order to confirm this, for a = 2, in Figure 6 we evaluated the expected ergodic rate 
of the system for different values of deployment radius r with Monte-Carlo simulation, where 
we used (13) for capacity. As we can see, the optimal radius r* = 360 is perfectly consistent 
with the optimal radius that the stochastic update method predicts. In general, as we can see 
in Figure 5 the lower bound approximation approach gives a better estimation of the optimal 
placement when we assume that the distributed ports are transmitting independently. 

C. Power Efficiency 

Optimal placement will increase the power efficiency of the networks. In this part we aim 
to quantify this amount for DAS(Q, 1), in comparison with the case where we have all of the 
antenna ports colocated in the center. In other words, for the case where we have the CSI only 
at the receiver and all antenna ports are co-located in the center of the cell, we set the per- 
antenna transmit power S so that the average SNR at the edge of the coverage region is 10 
dB. According to this transmit power for each value of a, we calculate the average ergodic rate 
C t (a) when all of the six antennas are co-located at the center. Then for each value of a, we 
solved Problem 4. 1 in order to find the optimal placement of the ports inside the coverage region 
and also the minimum power required to achieve the same average ergodic rate C t (a). Figure 
7(a) illustrates the ratio of the minimum powers required to satisfy C t in two cases. This gain is 
only due to the location optimization and as the figure indicates, the power gain increases with 
increased a. The same procedure can be repeated for the case where the CSI is available at the 
transmitter as well. Figure 7(a) also shows the gain that we get from optimally placing the ports 
when we have CSIT As we can see, the power gain is considerably smaller when channel side 
information is available the transmitter. In other words, optimally distributing the antenna ports 
in the case of having CSIT doesn't help much, whereas the gain is noticeable when channel 
side information is not available at the transmitter, but only at the receiver. Further power saving 
is possible by utilizing the single transmission strategy of Section V-B. Figure 7(b) shows the 
power gain we achieve from using the single transmission scheme in DAS(6, 1) in comparison 
with the case where all of the six antenna ports are co-located in the center of the cell and 
we only have CSIR. Note that in the case of having CSIT, co-located ports can use a coherent 
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transmission and reduce their transmit powers and yet achieve the same average ergodic rate. 
Therefore, the power gain from using single transmission strategy would be around 7.8 dB less, 
when the channel side information is available at the transmitter. 

D. Optimal layout is not always circular 

In this part we run the stochastic update method in order to show that placing the antenna 
ports in a circular layout is not always optimal. In particular we run the stochastic update method 
to find the optimal placement of DAS (6, 1) in an environment with path loss a = 5. Figure 
8(a) depicts the trajectory of the location of the ports through the iterations of stochastic update 
method, with random initial points. Yellow points in this figure are the initial locations we started 
with. As we can see the optimal deployment is not circular and one port should be located in 
the center. Figure 8(b) also shows the optimal layout for the case where we have N = 12 ports. 
As we can see, in this case, the optimal layout has no port in the center of the cell. 

E. Interference effects 

In this part we study the effect of interference in the optimal placement of antenna ports, as 
described in Section IV-C. For N — 7, L — 1, we observe that the optimal layout always has 
one port in the center and six other ports in a circle of radius r around it (see figure 1), where 
r depends on the interference coefficients 7, as well as the path-loss coefficient a. In this part 
we set all the interference coefficients from neighboring cells j — 1, 2, 6 to 77 = 7. We also 
let the stochastic gradient method find the optimal power allocation vector S* for antenna ports. 
Figure 9 plots the optimal radius r of the antenna ports around the central port, as a function 
of interference coefficient 7, for different values of path-loss a. As we see in the figure, the 
optimal layout shrinks towards the center of the cell as the interference coefficient 7 increases. 
For the power allocation, we observe that the optimal power allocation of all ports around the 
central one are the same. Figure 10 illustrates the ratio of the optimal power of the central port 
to the optimal power of the peripheral ports, as a function of a. 

F. Cell size 

Studying the effect of cell size on our antenna port placement optimization framework is also 
of interest. If we fix the transmit power of the antenna ports, scaling down the cell size by a 
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factor of K leads to an increase of the overall power consumption of the network proportional 
to K 2 . Therefore in order to have a fair comparison in this part, while scaling down the cell 
size, we also reduce the transmit powers by a factor of K 2 . On the other hand, if we neglect 
the near field effect of r 0 in (5), scaling down the distances reduces the path-loss power falloff 
proportional to K~ a . This leads to an increase of both the received powers from the antenna 
ports and interferers, by a factor of K a ~ 2 . Equivalently, we can say that scaling down the cell 
size by a factor of K reduces the noise variance a 2 in (14) proportional to K a ~ 2 . Hence for 
the path loss exponents a > 2 and for large values of K, we will have an interference-limited 
system. 

In order to illustrate the performance of optimally placed DAS as the cell size shrinks, we 
use the area spectral efficiency (ASE) metric in this part. This efficiency metric [17] is defined 
as the average ergodic rate of the user per unit bandwidth per unit area, i.e 

A(V , S) = C -^, (35) 

where P and S are placement and power allocation of DAS in a cell of radius R. Figure 11 
illustrates this metric for N = 3 single antenna ports, optimally placed in a cell with a varying 
radius 7 . The same figure also plots the area spectral efficiency of a system with N = 3 ports, 
randomly placed inside the cell. For both cases we scale the power proportional to the area of 
the cell in order to have the same power consumption per unit area. Also for both cases we 
assumed that the CSI is only available to the receiver. As we see in the figure, in comparison 
with random placement, the optimal placement of the ports in DAS has a significant effect on 
area spectral efficiency of the system when we have CSIR. Adding this result with the fact that 
even randomly placed DAS outperforms colocated antenna systems [8], we can conclude that 
optimally placed distributed antenna system can improve the spectral efficiency of the system 
significantly. 

VII. Conclusion 

We established the downlink capacity of a single antenna receiver for generalized distributed 
antenna systems where each port has multiple transmit antennas. We considered different as- 

7 Note that for N = 3 optimal power allocation is the same for all three ports 
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sumptions about the availability of the channel states at the transmitter together with a per- 
port power constraint. For the case where CSI is available both to the transmitter and the 
receiver, we showed that the optimal signaling for each of the antenna ports is beamforming, 
with the optimal beam weights proportional to the channel vector from the port to the receiver. 
For the case where the CSI is only available to the receiver but not to the transmitter, we 
showed that under a symmetric fading assumption, the optimal strategy for each port is to send 
independent signals with equal power from it's antennas. We have also presented a general 
framework for the placement optimization problem of distributed antennas in cellular networks. 
We concentrated on fully cooperative networks in which the base stations are connected together 
via an ideal backbone. We presented a fairly general framework with no implicit constraint on 
the antenna locations. As we showed in Section IV-C, our stochastic optimization framework 
is sufficiently general to incorporate interference as well as general performance metrics. The 
proposed algorithm also works for arbitrary coverage regions. Since we assumed that we only 
have one tier of interferers in our SINR calculations, an obvious extension would be to incorporate 
the impact of frequency reuse on optimal placement of the ports. Finally we should mention that 
our optimization framework can also be used in the optimization of relay placement in cellular 
networks. 

Appendix I 
Capacity with CSIR 

In this appendix, we first generalize the proof of [7, App. B] in order to show that the 
optimal covariance matrix Q* that solves (9) for a generalized DAS(N, L), should be diagonal. 
Then we prove that the power budget in each antenna port should be equally allocated among 
the antennas. We denote the submatrices of the Q, the correlation of the Gaussian transmitted 
signals of the mth and nth ports, as q mn = E [a;^a; n ] ixL 8 - First, note that tr (q nn ) must be equal 
to S n , otherwise, we can increase the diagonal elements of q nn until satisfying the constraint 
tr (q nn ) < S n , and hence increase the objective. The main issue is to find the correlation matrices 
q mn for m^n. 

8 We dropped the index 0 of the cell from the channel gain vector and transmitted signals 
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The objective function in optimization problem (9) can be expressed as 
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Now, the key observation is that by exchanging /ii by — /ii the optimization problem does not 
change due to the symmetric distribution of hi. Therefore, by flipping the sign of hi, the objective 
function does not change 
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where = {(to, n)\2 < m,n < N,m ^ n}, and Ki = {(m,n)\l < m,n < N,m ^ n}\Ki. 
Note that in (40) we flipped the sign of hi and since either to or n equals 1 in the set K\, we 
should flip the sign of the summation over K x . Using (39) and (40), we can write the objective 
function as 
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where equality in (42) holds iff -^^2 Kl h m q mn h^ = 0 for all h m ,h n , where(m,n) e -K'i. In 
other words, the equality holds if and only if q in — q ni — 0 for all n ^ 1. Repeating the same 
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procedure by flipping the sign for each of the channel vectors h m we can get the following 
inequality: 



log 2 



1 N 



H 

nn'i'n 



Qnnh 
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(43) 



where equality holds if q mn = 0 for all m ^ n. Similarly, for each of the antenna ports, because 
of the independence and symmetry of all NL channels in our model, we can repeat the same 
trick of flipping the sign of the channel gains for each antenna element to show that the off 
diagonal elements of the covariance matrix q nn should be zero, for all 1 < n < N. In other 
words we can rewrite (9) as 
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(44) 



where g nn (M) is the ith diagonal element of q nn . Now (44) is a convex optimization problem 



with the optimum solution of q n 



{1,1} 



S n /L (can be obtained from KKT conditions). One way 



to show that the power budget in each port should be allocated equally between the antennas 
is to check the KKT conditions for optimality of (44). The other way is to show that for any 
pair of antennas i ^ j in the n th port, replacing the power allocations q nn (i,i) and q nn (j,j) by 
their average power q nn = | {q nn {hi) + Qnn{j ■,])), the objective function in (44) increases. To 
see that, consider the simple case where N = I and we only have two antennas, i.e. L = 2. 
Assume that the optimal values of allocated power for the antenna ports that maximizes (44) 
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are q\ x {l, 1) and q\ x {2, 2). Now we can write the optimal value of (44) as 
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where A = ^(l, 1) — q, and in deriving (46) from (45), we used the fact that the channel 
statistics for both of the antennas in the port (h(l) and h(2)) are identical and independent. 
Therefore by swapping the index of antennas, the optimal value of the objective function should 
not change. The last inequality indicates that the optimal value of (44) is less than (47), which is 
a contradiction unless we assume that ^(l, 1) = q\ x {2,2), or equivalently A = 0. This proves 
that the power should be equally allocated to both antennas for the case of having a single port 
and L = 2. The proof can be easily generalized for any number of ports and also L > 2. 



Appendix II 
Capacity with CSIT 

Since logarithm is a monotonically increasing function of its argument, maximizing log(l + 
h Q h H ) in (11), is equivalent to maximizing the SNR = h Q h H . Therefore, the objective 
function that we want to maximize in this part can be written as 
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Now we want to prove that for all 1 < m < N and 1 < n < N, the following inequality holds 
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where for deriving the last inequality, we used the Jensen's inequality together with the covariance 
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the power constraint Y,i=i^[\\ x m(i) || 2 ] < Sm, we can conclude that J2i=i \\ h m(i)\\ \ZE[\\x m (i)\\ 2 ] < 
\\h m \\^/S^i, for all 1 < m < N. Replacing this in (49) results in 
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Then, from (48) we can write the following inequality for the objective 
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Now, note that the equality here holds if we replace the covariance matrix with 

Qlsit = [Qi,Q2, 1n] H [qi, 92, -, 9jv] , 



(51) 



(52) 



where q n = ^ h ^ y/S n , 1 < n < N. Since Q* csit also satisfies the per-port power constraint 
in (11), we conclude that Q* csit is indeed the solution of (11). 
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Fig. 1. Structure of distributed antenna system in a cellular setting 
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Fig. 2. Average ergodic rate of a user, randomly placed in the coverage region, as a function of 




Fig. 3. 
shapes. 



Three port placement within a hexagonal region of radius R. Location of the transmitter ports are denoted by star 




Fig. 4. Trajectory of the location of the ports in iterative algorithm with random initial points. Yellow points are the initial 
locations we started with 
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Fig. 5. Optimal placement radius r of design example 4.1 for lower bound approximation approach vs stochastic update method 
in hexagonal cell of radius R — 1000 




Fig. 6. Average ergodic rate of a single user in example 4.1 for deployment of antennas in different radius r 
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(a) Power gain without single transmission strategy 
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(b) Power gain with single transmission strategy 



Fig. 7. Power gain obtained from optimal placement of ports in design example 4.1 




(b) N = 12 

Fig. 8. Optimal placement of N ports in an environment with path loss a = 5. 




Fig. 9. Optimal radius r of the antenna ports as a function of 7 




Fig. 10. Optimal positioning of the ports vs. path-loss exponent for different interference coefficient 7 s 
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Fig. 11. Area spectral efficiency for a cell of radius R for path loss exponent a = 3, and interference coefficient 7 = 0.5 



